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We give a sequence of equivalent formulations of the ADE and ADE height models defined on a random 
triangulated surface: random surfaces immersed in Dynkin diagrams, chains of coupled random matrices. Coulomb 
gases, and multicomponent Bose and Fermi systems representing soliton r-functions. We also formulate a set of 
loop-space Feynman rules allowing to calculate easily the partition function on a random surface with arbitrary 
topology. The formalism allows to describe the critical phenomena on a random surface in a unified fashion and 
gives a new meaning to the ADE classification. 

Talk delivered at the Conference on recent developments in statistical mechanics and quantum held theory 

(10 - 12 April 1995), Trieste, Italy. 



1. INTRODUCTION 

During the past decade, stunning progress has 
been made in understanding and solving of sta- 
tistical models on two-dimensional random lat- 
tices. The interest in these systems arose mainly 
because of their interpretation as models of dis- 
cretized quantum gravity or, equivalently, "non- 
critical" string theories. The geometry of the lat- 
tice is considered here as an additional fluctuating 
variable. 

The simplest statistical model on a random lat- 
tice is the random lattice itself (the target space 
consists of a single point). This problem goes 
back to the problem of counting planar diagrams 
solved in the seminal paper |^ after having been 
reformulated in terms of a large N matrix inte- 
gral. More recently, some nontrivial statistical 
models have been solved on a random lattice fol- 
lowing the same principle: the Ising model J3| , the 
0{n) model the Q-state Potts model y], and 
the ADE height models jsj . Each such system is 
equivalent to an ensemble of coupled random ma- 
trices labeled by the points of the target space of 
the model. Besides the models listed above, there 
are many other solvable ensembles of matrix mod- 



*On leave from the Institute for Nuclear Research and Nu- 
clear Energy, 72 Boulevard Tsarigradsko Chaussee, 1784 
Sofia, Bulgaria 

tkostovOamoco . saclay . cea.fr 



els whose geometrical meaning still remains to be 
understood. 

The statistical systems on random lattice are 
interesting not only as toy models of discretized 
quantum gravity. In spite of their simplicity, the 
models on a random lattice still contain consider- 
able information about the critical behavior of the 
models on a regular lattice. Therefore they can 
be used as an heuristic instrument for studying 
complicated situations on regular lattices. For 
example, the dilute critical regime of the ADE 
models has been first discovered first on a ran- 
dom and then on a regular lattice 0. Fur- 
thermore, there are infinitely many multicritical 
regimes known on a random lattice, which are not 
yet considered on the regular lattice. We believe 
that the correspondence between the models on 
random and regular lattice merits to be under- 
stood on a deeper level and might help to estab- 
lish the missing connections between the various 
resolution methods used by now. 

In this talk I would like to review the con- 
struction and the solution of the ADE and ADE 
height models on a random lattice, as well as to 
present several equivalent formulations of these 
models as simple integrable systems. 

It is well known Q that the minimal two- 
dimensional conformal-invariant QFT are classi- 
fied by the simply laced Lie algebras (i.e., these 
of the classical series Ar, Dr, Eq, Et, Es). Each 
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of these theories can be constructed microscop- 
ically as a lattice statistical model ( a height 
model), whose local degrees of freedom (heights) 
are labeled by the points of the Dynkin diagram 
of the corresponding Lie algebra (|]. Similarly, 
the height models associated with the extended 
ADE Dynkin diagrams describe conformal invari- 
ant QFT with C — 1 and discrete spectrum of 
conformal dimensions. The height models in their 
dense (dilute) version can be mapped onto the 6- 
(19-)vertex model and solved using Bethe Anzatz 

(0); 

The height models can be defined on an ar- 
bitrary surface made of triangles or squares in 
such a way that the mapping onto a vertex model 
still holds. In the next section we will remind 
the definition of the height models on an arbi- 
trary triangulated surface. It requires supple- 
mentary Boltzmann factors associated with the 
points with nonzero local curvature on the sur- 
face. Then we describe briefly the loop gas rep- 
resentation and the loop diagram technique that 
follows. 

The models on a random lattice become topo- 
logical in the sense that their correlation func- 
tions do not depend any more on the distance. 
The only parameters left are the topology of the 
lattice, the volume and the length of the bound- 
aries and the boundary conditions. The local 
scaling operators can be introduced by shrinking 
a boundary with given boundary condition. 

In Sect. 3. the height models on random lat- 
tice are then reformulated as systems of coupled 
random matrices, or, in the eigenvalue represen- 
tation, as Coulomb gas systems. 

A complete set of observables in the height 
models on a random surface is given by the set 
of all loop amplitudes or, in terms of free bosons 
or fermions, the correlation functions of currents. 
A genus 5, n-loop amplitude is equal to the par- 
tition function of the model on a random surface 
with the topology of a sphere with g handles and 
n boundaries with given lengths. In Sect. 4. we 
present the loop diagram technique established 
in fisll and allowing to compute any loop ampli- 
tude on a surface with any topology. It will be 
formulated as a set of Feynman rules involving 
propagators, vertices (including tadpoles) of all 



topologies, and leg factors for the external bound- 
aries. A vertex of given topology factorizes into 
a fusion coefficient for the order parameters of 
the height model and an "intersection number" 
associated with the corresponding punctured sur- 
face. In Sect. 5. the partition functions are refor- 
mulated, using the vertex operator construction, 
as systems of free bosons or fermions (soliton r- 
fimctions). 

2. INTEGRABLE HEIGHT MODELS 
LABELED BY DYNKIN DIAGRAMS 

2.1. The target space 

Let X be the Dynkin graph a simply laced 
Lie algebra of rank r. Such a graph consists of 
a set of r nodes, labeled by an integer height 
X e {l,2,...,r}, and a number of bonds between 
nodes. Two nodes are called adjacent (~) on X 
if they are connected by a bond. The graph X is 
defined by its adjacency matrix 

j^xx' i if ^ ^ ^ , f-t\ 

\ 0, otherwise 

which is related to the Cartan matrix of the scalar 
products of the simple roots a^, by 

a^-a^' =2^:,. -A^^'. (2) 

Each Dynkin graph X can be extended to a 
graph X with ?' + 1 nodes by adding an extra 
node representing the lowest root. The graph 
X can be considered as the Dynkin graph of 
the corresponding affine Lie algebra. For exam- 
ple, the Dynkin diagram of the affine algebra 
[R = 1, 2, ...) represents a closed chain of i? -I- 1 
points. One can consider as well the case R — 
where the target space, which will be denoted by 
Aq, consists of one oriented link with identified 
ends. The connectivity matrix of this target space 
is Aoo = 2. 

In what follows we assume that the target space 
X is a Dynkin graph of an extended Dynkin 
graph. 

The eigenvectors 5^™^ of the adjacency matrix 
are labeled by the integer Coxeter exponents m 

^^x.'^M ^2cos^ (3) 

x' 
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where h is the dual Coxeter number of the Lie 
algebra. 

In order to simplify the notations we shall al- 
ways denote the eigenvector with largest eigen- 
value (the Perron- Frobenius vector) by Sx , 



^ _(si'\ A,D,E- 
^''~\S^^\ A,D,E. 

The normalized quantities 



X: 



(m) _ ^ 



(m) 



satisfy a closed (fusion) algebra 

(m) (m') ^ /7O „v(™") 

Ax A.X / J ■mm m A-x 

(m") 



(4) 



(5) 



(6) 



and orthogonality conditions of the form 

^x Xx Xx — Om.m' ■ \i ) 



E 



More generally, we define the genus-g fusion co- 
efficients 



Eo2-2g (mi) , , , (m„) 
^x Ax Ax 



(8) 



2.2. Definition of the height models 

Let us remind the definition of the height model 
with target space X on an arbitrary triangulated 
surface (two-dimensional simplicial complex) S 
j5|. Each height configuration on S represents 
a map S —>■ X such that the heights x{a),x{a') 
of any two adjacent sites a, a' of S are adjacent 
or equal in the target space. The partition func- 
tion of the height model defined on the surface S 
is equal to the sum over all height configurations 
S^X 



:f[s] =J2wis^ X) 



(9) 



with Dirichlet boundary conditions , i.e., con- 
stant height Xi along the connected components 
Ci of the boundary 95, imposed. The Boltzmann 
weight of each height configuration is a product 
of factors associated with the vertices a and the 
triangles A of 5 



Note that in the case of a regular triangular lat- 
tice, where exactly 6 triangles meet at each ver- 
tex, the Boltzmann weights associated with the 
vertices can be distributed among the adjacent 
triangles. In the case of a lattice with defects 
we have to take into account the local curva- 
ture and the the form (^0|) of the Boltzmann 
weights is more appropriate. The local Boltz- 
mann weights are expressed through the compo- 
nents of the Perron- Frobenius vector Sx as follows 



W,{x) = Sx 



W/s,{xi,X2,X'i) 



f\ A A 

2^2 3^3 ^X^iXi '^XiX2 



A 4 4 

'-'x 1X2-^X2X3-^X3X1 



r\ A A 

'-'X3X1 -^XiX2 -^X2X3 



(11) 



(12) 



2.3. The height models on a random sur- 
face 

The height model on a random lattice is the 
statistical ensemble of all height configurations 
S ^ X where the surface S is treated as as an ad- 
ditional dynamical variable. The partition func- 
tion of the model is defined as the sum over all 
surfaces S and all maps S ^ X satisfying certain 
conditions. Thus the height model on a random 
lattice can be also viewed as a problem of a ran- 
dom surface S immersed in the discrete target 
space X. 

Denote by Sf^'"^ ^ the ensemble of all trian- 
gulated surfaces S with genus 5, area A, and 
with boundary dS consisting of n loops Ci, C„ 
with lengths /i, Z„ For each surface S we define 
the partition function J-[S, xi, Xn] where Xi de- 
notes the height of loop d. The genus g, n-loop 
amplitude reads 



■g,A 



li.xi 



(13) 



The partition function of the height model with 
target space X is by definition the exponent 



Z[k,T, gkx 



(14) 



W{S -> X) = n W,[x{<^)] n W^a[x(A)]. (10) .F[«,r,5..J = E 



cres 



of the generating function for the loop amplitudes 
^ 1 1 



9=0 



A>0 n=0 
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m„ 



~/2 h,xi-. 



(15) 



The mapping of the height models onto a gas 
of loops yields a powerful combinatorial method 
allowing to write a complete set of equations for 
the loop amplitudes. Let us sketch the loop gas 
method. For given triangulated surface S, each 
height configuration defines a set of nonintersect- 
ing loops on the dual tri-coordinated graph. For 
example, if S is the regular triangular lattice, the 
loops live on the dual honeycomb lattice. 

The loops represent the domain walls separat- 
ing the domains of constant height x. The Boltz- 
mann weight of a loop configuration can be orga- 
nized as a product of factors associated with its 
connected domains and the boundaries of the sur- 
face ||] . A domain of height x contributes a factor 
S'^~^3~"', where n is the number of its boundaries 
and g is the enclosed genus. With the definition 
( p^ ) of the generating function, this rule applies 
also to the domains adjacent to the boundaries. 
The Boltzmann factors associated with a bound- 
ary of length I and height x is gix- 

The fact that the Boltzmann weights of the do- 
mains depend only on their topology allows to 
perform the sum over loops and the sum over 
surfaces in the opposite order. Namely, we first 
perform the sum over the shapes of the domains 
and afterwards the sum over the topologically in- 
equivalent loop configurations. The result of the 
summation can be expressed as a sum of con- 
nected graphs whose vertices and lines label cor- 
respondingly the domains and the domain walls. 
Let M^/^^"* ; be the number of of domains of genus 
g and n connected boundaries of lengths /i, Z„ 
(we assume that each boundary has a marked 
point). Then the partition function ( p^ is equal 
to the sum of all connected Feynman graphs com- 
posed from the propagator 



l + V l\ l'\ ' 



(16) 



and the tadpoles 



2g-l 



(g) 



W, 



ia) 



+ Sgfi gix 



(18) 



The large I asymptotics of the quantities 
is given by the loop amplitudes for the so 



(9) 



called "topological gravity" whose explicit form 
is known. The corresponding loop Feynman rules 
are derived in Section 5. 

The loop gas representation is encoded in the 
following Ward identities for the partition func- 
tion [| dH) 

n-2 



d 



dgn 



d 



k>3 



dgn~2-k,x dgkx 



E 



d 



d 



dgn-2-k,: 

d 



X 

d 



{k + k'y. 



1 



k\ k'\ T^+k' + l 



k,k'>0 

Z = {xeX,n= 1,2,...) (19) 



dgn-i+k,x dgk 

The coefficients of the genus expansion of ( [l9| ) 
give identities relating the loop amplitudes of dif- 
ferent genus. The leading term in the expansion 
in K leads to a closed equation for the disc ampli- 
tude Wx.i = dlogZ/dgk.xln^o 



Wn, a; = ^ kgk,xWn-2-k,x + ^ Wk,x 



k>3 



fe=0 



Wn-2-k.x+J2^""' E 



{k + k')l 



1 



kl k'l T'=+'='+i 



x' k,k'>0 

-l+k,x Wk',x] (a;eX,n = l,2,...).(20) 

3. ENSEMBLES OF COUPLED RAN- 
DOM MATRICES 

The crucial observation allowing to establish 
the equivalence with ensembles of random matri- 
ces is that the the vertices ( |T7| ) can be represented 
as the loop correlators in a Gaussian N x N ma- 
trix model with N = Sx/n 



the vertices 



yig) _ w(ff) A 

^x;h,...,l„ — '^*^ii,...,i„ I g 



2g-2+n 



,(n>2) (17) 



(9) 



g=0 

DMe-^'^'^- TrM'i...TrM' 



(21) 
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Therefore, if we associate with each height x a 
hermitian matrix variable Mx of size x 
where 



= — , 

K 



(22) 



the perturbative expansion of the matrix integral 



exp[Tr5:(-^M^ + ^g,.M^ 



1=1 



xex 



(23) 



around the gaussian measure will give the parti- 
tion function Z = (15) fl^ . One can check 
that the Ward identities ( [19D follow from the 
translational invariance of the integration mea- 
sure in (H^). 

The sum in the exponent with the cocfRcients 
( p^ is a simple logarithm and after shifting the 
matrix variables as 



^ + T/2 
we write the integral (|23| ) as 

Y[ I det(l (g) Mj; -I- M^' (g 



(24) 



-Try (M^) 



(25) 



where the coefficients in the expansion of the po- 
tential 



(26) 



In particular, the Aq model is dual to the 0{2) 
model on a random lattice and its partition 
function is given by the following N x N matrix 
integral 



[V,N] 



DM 



g~Try(M) 



det(l ig) M -I- M (g) i; 



(27) 



3.1. Coulomb gas representation 

The integrand in (|2^) depends only on the 
eigenvalues Xix,i = 1, 2, A^^^, of the matrix 
variable Mx- Therefore we can retain only the 
radial part of the integration measure dM.x ^ 
Hi '^^ix Th^{^ix - XjxY and write the partition 
function (g5|) as 

zx[v\N,-xex]= / nri'^^*.- 



Similar integrals have been introduced by the 
ITEP group as "conformal matrix models" [|l5|. 
To make the connection with the construction in 
fl^ j let us divide the nodes of the Dynkin diagram 
X into two groups (say, black and white) so that 
the nearest neighbors of each node have the oppo- 
site color and change the sign of the A-variables 
associated with the black points. Denote by 
the color (white (-I-) or black (— )) of the point 
X £ X. Then the integral can be written in the 
form 

Zx[V^.N,-x€X] = / \[ \{d\^ 
•' xexi=i 

g-y-(e.A„) -Q (A„-A,v)^"^-"^' (29) 

where we used the definition of the Cartan ma- 
trix (|). This integral can be interpreted as the 
partition function of TV = Nx Coulomb parti- 
cles with vector charges proportional to the sim- 
ple roots of the Lie algebra with Dynkin diagram 
X, restricted on a line and interacting with a com- 
mon potential. 

The integral over the eigenvalues (the positions 
of the charges) is well defined in the large N limit 
only if the potential keeps the white charges in the 
negative half-space A < and the black charges in 
the positive half-space A > 0. Then a nonsingular 
large A'' limit exists since the charges of the same 
color do not interact and the the black and white 
charges that attract each other are separated by a 
potential wall. The critical situation arises when 
the edges of the black and white charge densities 
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meet at the origin. The critical singularity is ex- 
plained in terms of the original statistical model 
by the dominance of triangulated surfaces with 
infinite area. 

In what follows we prefer to keep the repre- 
sentation ( ^ ) of the Coulomb system, which is 
more naturally related to our original problem, 
and which can be considered as a generalization 
of the partition fimction of the Aq model 



N 



(30) 



The existence of a nonsingular thermodynami- 
cal limit then implies that the distribution of the 
eigenvalues has its support along the negative real 
axis for all x. 

Therefore the integration over the eigenvalues 
Xix can be restricted from the very beginning to 
the interval A < 0. This restriction will mod- 
ify the nonperturbative sector but the genus ex- 
pansion of the partition function will remain the 
same. 

Note also that, since the point A = is out- 
side the support of the spectral density, one can 
consider a more general potential containing neg- 
ative powers of A as well. On the other hand, the 
algebra simplifies considerably if the even powers 
in the expansion of the potential V{X) are sup- 
pressed. 

It is sometimes more advantageous to consider 
instead of (EST) the canonical partition function 



A'x> 



Zx[V'',N.4. (31) 



We will assume in the following that the chemical 
potential //^ for the charges of type x is absorbed 
in the potential as a constant term. 

3.2. The partition function of the Ar model 
as a Predholm determinant 

Let us introduce the integration kernels 



^-l[V-(\) + V- (A')] 



A + A' 

(a; = 0, 1, x' ^ x + 1) 



(32) 



where the integration goes along the negative real 
axis. With the help of the Cauchy identity 



1 



det ■ 

y Ai + 



(33) 



ni<j(^»-^j)(M»-Mj) 
n.,,(A»+Mj) 

we can write the canonical partition function of 
the Ar model (whose target space X is a circle of 
R + 1 points) as a Fredholm determinant 

Z^^ = Detil + K01K12K23-KR0). (34) 

In particular, the canonical partition function 
for the Aq model is given by 

Z, dK det K{\, \, ) 



7v=o ■ i=l 



Dct(l + X) 



(35) 



where 



i^(A,A') = 



-i[y(A) + y(A')] 



A-^ A' 



V{\) 



2n+l 



(36) 



Fredholm determinants like ( |35| ) give the cor- 
relation functions in various statistical or QFT 
models. The simplest case is the two-spin cor- 
relation function in the Ising model satisfy- 
ing the Painleve III equation. More recently, this 
Fredholm determinant has been considered in the 
context of self-avoiding polymers on a cylinder 
II. 

The representation (|34| ) of the partition func- 
tions of the Ar models leads to the problem of 
diagonalization of the kernel (|3^). This problem 
has been solved explicitly only in the case of a 
potential V{z) = iiA + t-iA"^ An original 
method to find the eigenfunctions in this case has 
been suggested by M. Staudacher |2^. Unfortu- 
nately the method does not work more compli- 
cated potentials. The idea is to look for a differ- 
ential operator commuting with the integral op- 
erator represented by the kernel K and therefore 
having the same spectrum of eigenstates. If we 
introduce the parametruzation A = — e" (—00 < 
u < 00), then the kernel (^) becomes 
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with potential 

= V^(e") = -^-<ie" -f_ie-" (38) 

and nonrestricted homogeneous measure du. It is 
easy to check that the hnear operator with kernel 
(^7|) commutes with the differential operator 



where 



H 



v'{u)f 



1 

and therefore has the same set of eigenstates 



(39) 



3.3. Loop equations the ADE and ADE 
matrix models 

The resolvent 



WJz) = lim (V 



N^oo '■^ — ' Z — Af 
i—1 



(40) 



represents a meromorphic function of z with a cut 
along the support [a, 6] (a < 6 < 0) of the classical 
eigenvalue density and behaves at infinity as 



1 



(41) 



The loop equations ( |l^ ) can be derived in the 
matrix model from the translational invariance 
of the matrix integration measure. The classical 
{k —>■ 0) loop equation for the resolvent Wx{z) 
reads [|j 

Jc_ ^TTl Z — W 

+ ^A^^V:,/(-u;)] = (42) 

x' 

where the contour C_ encloses the point z and 
the cut [a, b] of Wx{z) and leaves outside the cut 
[-b,-a] of VF,(-z). 

The exact loop equations are obtained from 
( ^2[ ) by replacing Wx{z) with the loop insertion 
operator 



dV^{z) 



(43) 



and represent the following quadratic differential 
constraints on the partition function 



/ dw 1 
1(2 27ri z — w 



T^''\w) Zx = 0, xeX (44) 



r(^)(z) 



dV'izy 



W(z) 



dV'{~z)J dV^{z) 



Let us introduce the current-like field 



dVx{z). 



(45) 



(46) 



dV^{z) 

where the covariant components Vx{z) of the po- 
tential are defined by 

y^(z) = 2F,(z)-^A-^'l4,(-z). (47) 



Then the operators ( |45| ) can be written in the 
form 

r(-)(z) = a$,(z)[a$,(z) -f ^ A--'9l.,,(-z)]. (48) 

x' 

The classical expectation value d^x{z) of the 
current (|4^) is related to the resolvent (^ ) by 
d^x{z) = Wx{z) — dVx{z) and is completely de- 
termined by the conditions 

$,(A + ^0) + $,(A - lO) - ^ A""'$,.(-^) 

= 0, (a < A < 6); (49) 

*.(^) = ~Vx{z) +Nx\nz + 0(i), (z ^ oo), (50) 

which are obtained by taking the imaginary part 
of the classical loop equation (|4^). The first 
of these equations means that the meromorphic 
function 



$^(z) = 2$,(z)-^A^-'4,,(-z). 



(51) 



has vanishing real part along the cut [a,b]. 

The constraints (|4j) form a set of mutually 
commuting Virasoro algebras, one for each point 
X of the target space. One can interpret the Vi- 
rasoro constraints for given point x as the loop 
equations for the topological gravity with 
potential depending on the fields associated with 
the adjacent points x' . In this way the solution of 
the topological gravity appears as the mean field 
problem in our models. In the next section we 
formulate a collective field approach utilising the 
solution of the topological gravity. 
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4. LOOP DIAGRAM TECHNIQUE 

The tree (genus zero) amplitudes can be cal- 
culated with the effective collective field theory 
in which the eigenvalues of the random matri- 
ces are considered as a continuous classical liq- 
uid. The quantum fluctuations of the liquid in- 
teract in a complicated way but the interaction 
potential can be determined exactly. The strat- 
egy is based on the exact solution of the mean 
field problem, namely the integral over the posi- 
tions \ix {i = 1, Nx) of the charges of type x 
for fixed distribution of the other charges. This 
leads us to a one-matrix model for generic po- 
tential whose solution is known. The free energy 
of this effective one-matrix model gives the inter- 
actions of the vibration modes of the eigenvalue 
liquid. 

Let us introduce collective field ^x{z) and the 
corresponding lagrange multiplier field V^{z) by 
inserting the identity 



1- J V^VVYl exp(^^ 



-F,(z)-^ln(z- A„)]y"(z)) (52) 



1=1 



in the r.h.s. of (Eq). The integration with respect 



to the Xi 



1, ...Nx, yields the mean field free 



energy defined by 



|nrfA,e-^(^')n(A.-A,) 

2—1 2<J 



(53) 



After this change of variables the original parti- 
tion function becomes a functional integral over 
the collective fields $ and V with nonrestricted 
homogeneous measure 



Z = y V^VV e-'5[*-^"'^-e^l, (54) 



2TTi 



27ri 



ln(z + z'). (55) 



The dependence on the potential V^{z) is 
through the asymptotics of the field $j; at z —> cx3. 



The string propagator and vertices are ob- 
tained by expanding the effective action (|55| ) 
around the mean field <i>^, determined by the 
large N saddle point equations. The dependence 
on the string coupling constant k — Sx/Nx is 
through the genus expansion of the effective po- 
tential F[T>^] = J2g>oNx~^^T^^'>[V% It is con- 
venient to make one more change of variables and 
introduce the field ^x{z) defined by 



d^x{z)^- 



dV^{z) 



-dV^z) 



(56) 



as an independent fluctuating variable. The fields 
^x{z) and ^x{z) have the same classical values 
^xiz) but on the quantum level they play differ- 
ent roles : the external source is coupled only to 
the field ^xiz) while the interactions affect only 
of the field <hx{z). 

It is of course possible to consider the disconti- 
nuities 

(bx{X) = --^^xW, 0.(A) = --3$,(A) (57) 

as independent functional variables. In any case, 
the functional integration measure in (55) is de- 
fined only after introducing a mode expansion for 
the two fields. 

4.1. Saddle point 

Upon an appropriate rescaling of the A- 
variable, the support [a, b] of the classical spectral 
density p^(A) = — id3<&^(A)/dA can be taken to 
be the interval [— L, —1] on the negative real axis. 

The solution of the saddle point equation in the 
continuum limit L — s- cx) reads, up to an arbitrary 
normalization, 



X = ADE 



dz 



Kiz) 



Sx (z + + {z - yG^~^Y 



27r/3|sin7r/3| 



where 



/3 = 1 ± - + 2m 
h 



(58) 



(59) 



and the string coupling constant scales as k ~ 
^i+/3^ The different branches (±,m) correspond 
the different critical regimes of the height model 
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(for details see [g| and Q). The regimes (—,0) 
and (+, 0) are sometimes referred as dense and 
dilute phase of the height model. The central 
charge of the corresponding conformal field the- 
ory is 



C = 1-6 



/3 



(60) 



The saddle-point solution for the case X — ADE 
is obtained as the limit /3 ^ 1 of (|58|). 

Our aim now is to define the Hilbert space of 
one-loop states, choose a complete orthonormal- 
ized set of eigenstates of the quadratic action, 
and finally express the interactions in terms of 
the mode expansion of the collective fields. 

4.2. Gaussian fluctuations 

It is consistent with the perturbative expansion 
to assume that the fluctuating fields are again 
supported by a semi-infinite interval, but its right 
end can be slightly displaced with respect to its 
saddle point value —1 due to fluctuations. 

Inserting the genus-zero contribution to the 
mean field free energy ( ^3|) 

f(")[\4] = -iy5R$,d3l>, (61) 

we find for the tree-level action 

= / {3?*.(A)d3l>.(A) 

+ [^ A--'$^,(_A) - 23fi$,(A)]d3$,(A)}. (62) 



The tree-level effective action (|62|) can be split 
into gaussian and interacting parts 



= S 



free 



(63) 



n>3 



where iS^*^'"-* is the genus-zero n-loop interaction. 

The gaussian fluctuations of the collective fields 
are those that do not shift the edge of the eigen- 
value interval; the fluctuations displacing of the 
edge are taken into account by the nongaussian 
terms in the effective action that represent the 
tree level n-loop interactions. 

We therefore define the Hilbert space Ti. as the 
space of real functions defined on X x [— cx),— 1] 



or, equivalently, the space of analytic fields having 
a cut along the interval [— cx), — 1]. 

For the purpose of diagonalizing the quadratic 
action it is very useful to introduce the map 



z(r) = coshr 



(64) 



transforming the space of meromorphic functions 
in the z-plane cut along the interval [— oo,— 1] 
into the space of entire even analytic functions of 
T. In the following we will keep the same letters 
for the analytic fields considered as functions of r 
and denote $(r) = $(z(r)). The disc amplitude 
dsa) as a function of r reads 



dz 



Sx cosh (3t 



K 27r/3| sin7r/3| 



(65) 



The operations 3 and 3? in the z-space become 
finite-difference operators in the r-space. Denot- 
ing d = 9/9r),we can write 



3$(-2;) sin7r9<i>(T), 
m{-z) = cos7r(9$(r) 

It is clear that the plane waves 

{x,T\m,E) ^Si'^^e'^^ 



(66) 



(67) 



form a complete set of (delta-function) normal- 
ized wave functions diagonalizing the quadratic 
action The latter is given by the expression ( |6^ ) 
where the integration is restricted to the interval 
A < 1. According to (|66|), the Fourier components 
of the fields tj) and (f> are related to these of <t and 
$ by ^^^){E) = isinh(7r£;)|.(„)(£;), ^^^){E) = 
i sinh(7ri?)<I)(^)(i?). The action (|6|) reads, in 
terms of these fields. 



5" 



p—m/h 



El arj 
L 2^ 



dE r t:E cos np 



sinh t:E 



-(<^-20) 



ttE cosh ttE 
sinh irE 



(68) 



By inverting the quadratic form in (|6^) we find 
the propagators in the {E,p = ^) space 



E,p) = G'f'HE,p)^GiE,p), 



G't"t'{E,p)^GiE,p)-GiE,-) 



(69) 
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where 



sinh t:E 



1 



nE cosh ttE — cos irp 
2 ^ 1 



7r2 ^ ^2 ^ (p + 2n)2 



(70) 



n— — oo 



The three propagators have the following dia- 
grammatic meaning: 

q4>4> is associated 

with the external (internal) lines of a Feynman 
diagram, and G"^*^ is the genus-zero loop-loop cor- 
relator. 

4.3. Interactions 

In order to compute the interaction potential 
we need the explicit answer for the one-matrix 
integral. 

The free energy F[V] of the one-matrix model 
(^) is most easily expressed not in terms of the 
potential V{z), but in terms of the field 9$(z) = 
dF^°ydV(z)-^dV{z) (111 - H). Themeromor- 
phic function d^{z) has vanishing real part along 
its cut [a = —oo,b] and therefore can be expanded 
as a series in positive half-integer powers of z — b. 
(Note that b = b[V] depends on the potential!) If 
zo is a point close to b, the field can be neverthe- 
less expanded as a series in the positive and nega- 
tive half-integer powers of z — za- The free energy 
F can be expressed in terms of the coefficients of 
the nonsingular at z = zq part of the expansion. 
The coefficients are given (up to numerical fac- 
tors) by the linear functionals {k\^)zo,k = 0, 1, ... 
defined by the generating function 



dz d^'^{z)/dz 



2ni ^ z — zq — 2u 



(71) 



OO 

fc=0 

Sometimes they are called KdV coordinates of the 
field <&. The KdV coordinates at different expan- 
sion points are related by 

k 

(72) 



(*|n).o+2„ = ^($|n + /c),„^ 

fc>0 



The choice of the expansion point zq is a mat- 
ter of convenience. The most compact expression 
corresponds to the choice zq = b, where b is the 
position of the branchpoint of ^{z), imposed by 
the condition 



In this case the explicit expression for of F[^] 
reads il 



F[<i>]=F(«)[$]-lln(l|<i>) 



2-2g 



E 



25-2+n>0 



{^1 ■ ■ ■ ^n}c 



ki,...,kn>2 

fclH \'kn=3g-3 + n 



(fci|$)b (fc„|$)j 



(74) 



The coefficients {fci • • • fc„}g {ki + ■ ■ ■ + kn = 
Sg — 3-|-n) are the the intersection numbers on the 
moduli space A4g,n of algebraic curves of genus g 
with n marked points | |23[ | . The intersection num- 
bers can be obtained from a system of recurrence 
relations equivalent to the loop equations ||l[ . In 
particular, the genus-zero the intersection num- 
bers coincide with the multinomial coefficients 



{^1 ■ ■ • fcn}o 



(fci 



fci! •••/=„! ' 

/Cl + • • • + fc„ = 71 — 3. 



(75) 



The genus g term in the expansion of the free 
energy is given by the restricted sum ([tJ) with 
ki + ■ ■ ■ + kn = 3g — 3 + n, which contains only 
finite number of terms. 

The expression of the free energy as a func- 
tion of the KdV coordinates {n\^)zg with fixed 
zq has the same form but the sum over ki, ...Tkn 
should be taken in the range ki > since the 
coordinate (0|<I>)zo no longer vanishes. Note that 
F is singular at <I> = and has to be expanded 
around some nontrivial background ^'^{z). The 
simplest nontrivial background to expand about 
is the "topological" background 



1 2 
V23 



,3/2 



(76) 



with only one nonvanishing coordinate (^'^ll)^;, = 
1/k. a generic potential can be considered as 
a deformation about this topological point. For 
deformation parameters tk — tk{^) defined as 



(fc|$),„ = -{Sks~tk) 



(77) 



($|0)6 = 0. 



(73) 



the genus-g term of the free energy reads, in terms 
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of tk, 



ki kn>0 

felH |-fc„=3g-3 + Ti 



n>0 

{fcl • • -fcnlg tfej • • -tfc^. (78) 



Now we are ready to find the vertices (|T^) for 
the genus g, n-loop interactions {2g + n—2>0). 
We have to expand the effective potential 
around the background ( ^8|) and taking zq = —1. 
Alternatively, we can expand around the topolog- 
ical background (|7q ) with zq ~ —1 and consider 
the background ( |58| ) as a perturbation. We prefer 
to do the latter because in this case the interac- 
tion vertices do not depend on the background. 
Then the deformation parameters tkx — tk{4'x) 
are defined by 



<5fe,l — tkx — 



(79) 



(Once the expansion point is — —1 is fixed, 
we will write (fc|*) instead of (A;|*)_i.) We know 
the explicit form of the propagator and the tad- 
pole in the base of plane waves, while the interac- 
tion is expressed in terms of the KdV coordinates. 
Therefore we need the transformation matrix be- 
tween these two mode expansions. The generat- 
ing function for the KdV coordinates is given, in 
the T space, by 

1-°° drC0STTdr4>(T) 



Y^-ikm=l dr 



k=0 



(80) 



cosh 



For (/)(t) — sin Et the integral is a Legendre func- 
tion and its expansion in u gives for the KdV co- 
ordinates of a plane wave 

{k\E)^ nEUkiiE), /c = 0,1,2,... (81) 
where 

Uo{a) = 1; 

^'^i'^) = i^Elij + h'-a% fc-1,2,... (82) 



3=0 



Eq. (^Tj) means that the linear functional (fc| acts 
in the space of odd functions (^(t) — —(f){—T) 
smooth at T = , as 



The KdV coordinates of the string background 
([5^ ) are obtained using the identity {k\d(j)/dz) — 
\ (fc + l|0) 

(fc|,^S) = -^ n,_i(/3), fc = 2,3,.... (84) 

The interaction vertices for the field are ob- 
tained by expanding the free energy of the one- 
matrix model around the deformation parameters 

characterizing the saddle-point solution. From 
(|8|) we find 

tl=t\= 0, tl = ~Ilk-i{P), k = 2,3, ...) (85) 

After going to the momentum space the factor 
5*^^"^+" in the genus expansion of yields 
the genus g fusion coefficient (||) and the vertex 
factorises into a product of an intersection num- 
ber and a fusion coefficient. The KdV represen- 
tation of the internal propagator is given by the 
symmetric matrix 

Gkk'{p)^mk'\G^^m; (86) 
Goo(p) = -ni(l-p), 

Goiip) = Gio{p) = -Tl2{l-p), 
Gn{p) = -202(1 -p)- 203(1 - p), ... (87) 

Finally, the external line factors read 

sinh(l — p)t 



{k,p\G^AT,p)=Ilk{dr) 



sin T:p sinh r 



To summarize, we have found the following 
loop Feynman rules: 
Propagator : 



Gkk'ip), P=-f^ 
Tadpole : 

T/(0) - ^fcr 

Vertices : {g + 2n > 2) 

y{9) 
ki ,rni , . . .,kn jTUn 

— ^^^{ki • ■ ■ kn}g C^l __^-^^. 



(89) 



(90) 



(91) 



(fc|l>) =^[nfe(aO dr4>{r)]r=0- 



(83) 



Let us give several examples. 
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(i) The three-loop genus-zero amplitude 

($,„i(zi)$„,(z2)«'m3(2:3)>0 = 



771177127713 



n 



sinh(^)r, 
s— 1 



(92) 



(m) The genus-one one-loop amplitude 

($,„(z))i - ^kUi}ini(a.) + ci„{02}it§ 



sinh(^i^)T 

^ 771 ^ 

sinTT^ sinhr 

h 

1 /m' 



sinh(V)r 
2 y J sin ^ sinh r 



(93) 



(im) The four-loop genus-zero amplitude 



9t2 



E 



/m 1\ /m 3 
V7i^ ~ 2/ vX ~ 2 



(-<0 ^0 



A sinh(l-^)r, 
sinTT^r sinhTc 



(94) 



The expression for the genus 1 tadpole (|9^) is in 
accord with the recent calculation by B. Eynard 
and C. Kristjansen in the 0{n) model. (In this 
case one has to put ■nmjh = arccos(n/2). ) Eq. 
( |9^ ) reproduces in the limit oo the four-loop 
amplitude found in ref. |24]. 



an established fact. It has been shown by M. 
Douglas that the chains of random matrices 
with nearest neighbor interaction (TrMxMx+i) 
are described by the A series of generalized KdV 
flows in the classification of Drinfeld and Sokolov 
p6|| . Moreover, as it was noticed by Di Francesco 
and Kutasov |2^ , a particular case of the D series 
reproduces the known qualitative features of the 
D models of matter fields coupled to gravity . 

It remains an open problem, however, how to 
construct the hierarchies associated with the ra- 
tional conformal theories coupled to gravity in a 
unified fashion, including the exceptional ones. 
We believe that the ADE and ADE matrix mod- 
els can be used as a starting point to solve this 
problem. The logic is the following. Each matrix 
model partition function in the form of Coulomb 
gas defines a bosonic vertex operator construc- 
tion. Then the associated hierarchy of soliton 
equations can be written as a system of Hirota 
bilinear equations |30| following the general pre- 
scription proposed in In the case Aq this 
is KdV hierarchy (just as in the case of the Ai 
model, but of course with different boundary con- 
ditions). 

Below we show how to make the first step, 
namely the bosonic vertex operator construction. 
We restrict ourselves to the A and A series, but 
the generalization to the other cases is evident. 
Then we give the fermionic representation for the 
A and A series. The boson-fermion correspon- 
dence allows to establish a fermionic representa- 
tion for the A, A, D and D series, but not for the 
exceptional [E and E) cases. There should be 
certainly a connection between the hierarchies as- 
sociated with the A and D Coulomb systems and 
the Lax representations by Drinfeld and Sokolov 
for the A and D series IM. 



5. VERTEX OPERATOR CONSTRUC- 
TION 

The relation between the Coulomb gas systems 
and the integrable hierarchies of differential soli- 
ton equations gives another powerful approach to 
study these systems, especially in application to 
nonperturbative phenomena. 

The relevance of the integrable hierarchies to 
the critical phenomena on random surfaces is now 



5.1. The case Aq 

Before presenting the general construction let 
us consider in detail the simplest nontrivial exam- 
ple, namely the Aq model whose the target space 
is a circle with one point. 

Let us first note that if the integration with re- 
spect to A is replaced by contour integration with 
respect to the complex variables z.^, the partition 
function will still satisfy the same loop equations 
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(ii- 

We start with the construction of the bosonic 
field representation of the canonical partition 
function 



N 



-V(z,) '■'■i<3 

n.,(- 



(95) 

I I . . ( — z.- — Z. I 

N=0 4=1 

Introduce the bosonic field ip{z) with mode ex- 
pansion 



(f(z)~q+p\nz + 'S^—z " 



n#0 



[Jn, Jm] = nS,n+n,o; [P, 4] = 1- 

and the vacuum state |^) defined by 
J„|0=0, (n>0); p\l)^m- 



(96) 
(97) 
(98) 



The associated normal ordering is defined by 
putting J„,n > to the right. Following ref. 
31 1 we define the antisymmetric field 



$(z) = ^(z) - (^(-z). 



(99) 



The field $(z;) and the vertex operator : e*'^-* 
satisfy the OPE 



$(z)$(z') =: $(z)$(z') : 
+21n(z-z') -21n(z + z'), 

. g*(2:) .. g*(2:') _ ~ 

(z + z')2 
Define the Hamiltonian 

n 

and the operator 



(100) 



Gao = exp (^e'' (j) 



dz 
27ri 



(102) 



(103) 



It is easy to see that the vacuum expectation 
value 



To[t] = (0|e^[*lG^jO) 



(104) 



is equal to the partition function ( |95D with poten- 
tial 



V{z) 



2n+l 



(105) 



The fermionic representation of the parti- 
tion function of the Aq model follows from the 
bosonization formulas 



V'(z) 



d(p{z) —: tp*{z)'ip{z) 
where the fermion operators 



(106) 



(107) 



satisfy the modes in the expansion of the anti- 
commutation relations 



(108) 



The operators ( |102| ) and (103) are represented by 



n>0 



G = exp 



=^ ^ :V.(z)^*(-z) 
27ri 



(110) 



and the vacuum states with given electric charge 
/ satisfy 



= (^IV^; = (r>0 



(Ill) 
(112) 



Following the line of arguments of ref. [£8|, 
we can consider the canonical partition function 
( p^ ) as the large iV limit of iV-soliton solutions of 
the KdV hierarchy of integrable differential equa- 
tions (where t\ is the space variable and ^3,^:5, .. 
the time variables). Therefore the partition func- 
tion (|9^) represents itself a r- function of the KdV 
hierarchy. 

More generally, one can define 



rat] = (?|e^WGAolO- 



(113) 



for any integer I but in fact there are only two dif- 
ferent r-functions corresponding to I = 0{mod 2) 
and / = l(mod 2). In terms of Fredholm determi- 
nants 



Tl = dct(l + {-YK) 



(114) 
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where the kernel K is defined by (|3q). 

The KdV and mKdV hierarchies of differen- 
tial equations are generated by the Hirota bilinear 
equations 

n 

Ti(i„--z-")T,,(C + -z-") = 0. (115) 
n n 

Namely, the second derivative 
dtj 

satisfies the KdV hierarchy of differential equa- 
tions, the first of which is the classical KdV equa- 
tion 



" 2^ log TO 



(116) 



du du 

Further, the function 

d . To 
v = — log — 
dti Ti 

satisfies the modified KdV equation 



(117) 



(118) 



dv 



dv 



(119) 



and u and v are related by the Miura transforma- 



tion 



dv 



(120) 



Finally, the bosonic and fermionic representa- 
tions generalize trivially if we add negative odd 
powers to the potential, 



ri[t] 



10, (121) 



and the r-functions (|l2lj) solve the affine sinh- 
Gordon hierarchy (with t_i,<_3, ... as time vari- 
ables). The lowest order equations are 



^ — — - sinh 26, 



92 



dtidt. 



logri = 



.20 _1 



(122) 
(123) 



(0 = log^). 
The relation between the Fredholm determinants 



(114) and the mKdV and sine-Gordon hierarchy 
has been established directly in | |T9| ]. 



5.2. Bosonic construction in the general 
case 

The idea is to construct the operators 
<^^{z), X € X, associated with the roots of the 
corresponding classical Lie algebra, with OPE 

$^(z)$^'(z') =: $^(z)$^'(z') : 
+2^^/ ln(z - z') - A^^' ln(z -|- z'). (124) 

We will restrict ourselves to the cases Ar and Ar ; 
the other simply laced algebras can be considered 
in a similar way. 

We will use the fact that the simple roots 
[x = l,...,r) and the lowest root cS° of A,. 
can be represented as differences of orthonormal 
vectors ei, ...,6*^+1 in H'^^ 



Cx-e^+ii (a; = 0, ...,r) 



a =6^+1-61. 



(125) 



With each vector Ca (a — 1, 2, r + 1) we asso- 
ciate the bosonic field i^aiz) having mode ex- 
pansion 



Va{z) = (7a + Pa In 2 -I- 

with 



n 



\Jan:Ja'm\ — ^^^a^a' ^n-\-7n,Qi 



[Pa,qa' 



(126) 

(127) 
(128) 



The bosonic Fock space JF^i ^ = {h, ■■■Jr+i} is 
generated by the action of the negative modes on 
the vacuum state \l) — e^^'^^'^\0) satisfying 



Pa 10 =^a|0; ^a«|0 =0, n > 0. 



(129) 



The state (/| is defined similarly, with the normal- 
ization {l\l') = Y[a^la,l'a- 

The fields 

<i>'={z)^(p^{z)-(p^+i{-z) (a; = l,...,r) 

$"(z) = (^,+i(z)-^i(-z) (130) 



satisfy (124) and therefore the corresponding ver- 
tex operators have OPE 

(z - z'f'^^' 



{z + z')^-.-'+^+^- 



g<I.-(z)g*=(z') 



(131) 
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Define the Hamiltonian 

r+l oo 

] ^ ^ ^ ^ ^ ianJan- 
a—1 71—0 



(132) 



and the operator creating the electric charges 
with fugacity e'' 



G; 



(133) 



Then the vacuum expectation value 
Tg[t;A.]^(0|e«[*lG|0) (134) 



is equal to the canonical partition function (28) 
for the model Ar with potential 

V^z) = -Ai-5]z"[i.„- (-)"Wi,n]. (135) 



More generally, one can define the r-functions 
associated with vacuum states with different vec- 
tor charges I and I' such that li + ... + Ir+i = 

+ ... + /;+! 



(136) 



The bosonic representation of the model is 
constructed in the same way. The canonical par- 
tition function is given by 



where 



Hit 



Ga^ 
and 



^ / 977-,- 



exp 



27ri 



(137) 



(138) 



QX 

, a; 

I^+l.O = 1^ — — 
K 



5^ = 2S':,-^:,:,/S':,' =4sin^ 



/ TT 



\2h 



(139) 



5.3. Fermionic representation of the Ar 
and Ar models with r > 2 

The fermionic representation of the partition 
functions follows from the bosonization formulas 



(140) 
(141) 
(142) 



Let us introduce the fermion operators 
i'r,a,ipr,a (r € Z 5, a = 1, R+l) Satisfying 
the anticommutation relations 

[lpr,a,1p*s,b\+ = SrsSab- (143) 

Define the zero-particle states | — oo), (— oo| by 



1pr,a\ - 00) = 0, (-CX)|V^*^ = 

and the states \l) = \li, ...,lr+i), 
{h,...,lr+i\, (ZaeZ)by 

io=nncai-oo) 

a r<la 

(/i=(-ooinn^-.- 

a r<la 

These states satisfy 

m.r.a = iWr^a =0 (r > Z,) 
V'r.alO = CalO =0 (r > L). 

Define the current operators Jn.a 

Jn,a = • '^r-n,aA,a ■ (n G Z) 

r 

where 

: V'*,q'0s,& := V'r.aV'sb - {0\llJr,a'^sb\0), 

or, equivalently, 

Ja(z) = -^"'^ =■ ^a(^)^a(z) 

where 

'/'a(2:) = X! A,aZ~''~K 



(144) 

(^1 - 



(145) 

(146) 
(147) 

(148) 

(149) 
(150) 



(151) 



Introduce the Hamiltonian H by (132) such that 



'>Pa{z)e' 



''Paiz) 



raiz)e- 



-H{t) 



rE:..*"."-"V':(^)- (152) 



Then the operators (|l33|) an d (138) in the defini- 
tion of the T functions (134) and (137) are repre- 
sented by 



Ga^ =expe^y : Vx(z)^*+i(- 



-z) 



+ :^r+i{z)ri{-z): (153) 
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and 



dz 

Zm ^ — ' 

x = l 



.(154) 



The T-functions for the Ar and Ar models 
satisfy the Hirota bihnear equations ||3^] foUow- 
ing from the commutatativity of the operators 
S = J2af dzipaiz) <E> i^liz) and G ® G. Bihnear 
equations for the D and E series can be written 
as weh but they have more comphcated form |29| . 
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